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Abstract —In this paper we propose distributed dynamic 
controllers for sharing both frequency containment and restora¬ 
tion reserves of asynchronous AC systems connected through 
a multi-terminal HVDC (MTDC) grid. The communication 
structure of the controller is distributed in the sense that only 
local and neighboring state information is needed, rather than 
the complete state. We derive sufficient stability conditions, 
which guarantee that the AC frequencies converge to the 
nominal frequency. Simultaneously, a global quadratic power 
generation cost function is minimized. The proposed controller 
also regulates the voltages of the MTDC grid, asymptotically 
minimizing a quadratic cost function of the deviations from 
the nominal DC voltages. The results are valid for distributed 
cable models of the HVDC grid (e.g. 7r-links), as well as AC 
systems of arbitrary number of synchronous machines, each 
modeled by the swing equation. We also propose a decentralized, 
communication-free version of the controller. The proposed 
controllers are tested on a high-order dynamic model of a 
power system consisting of asynchronous AC grids, modelled 
as IEEE 14 bus networks, connected through a six-terminal 
HVDC grid. The performance of the controller is successfully 
evaluated through simulation. 

I. Introduction 

Power transmission over long distances with low losses 
is an important challenge as the distances between genera¬ 
tion and consumption increase. As the share of fluctuating 
renewables rises, so does the need to balance generation 
and consumption mismatches, often over large geographical 
areas, for which high-voltage direct current (HVDC) power 
transmission is a commonly used technology. In addition 
to offering lower cost solutions for longer overhead lines 
and cable transmission a the controllability of the HVDC 
converters offers flexibility and means to mitigate problems 
due to power fluctuations from renewables. Increased use 
of HVDC technologies for electrical power transmission 
suggests that future HVDC transmission systems are likely to 
consist of multiple terminals connected by HVDC transmis¬ 
sion lines, so called multi-terminal HVDC (MTDC) systems 

El. 

The fast operation of the DC converters enables frequency 
regulation of one of the AC grids connected to the HVDC 
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link. One such example is the frequency regulation of the 
island of Gotland in Sweden, which is connected to the 
much stronger Nordic grid through an HVDC cable a. 
By connecting multiple AC grids by an MTDC system, 
enables frequency regulation of one or more of the AC grids 
connected. Traditional AC frequency controllers and HVDC 
voltage controllers do however not take advantage which the 
increased connectivity of the grids brings. Rather than sharing 
control reserves, each AC area is responsible for maintaining 
its own frequency in an acceptable range 0, which reduces 
the need for frequency regulation reserves in the individual 
AC systems 0 , & A challenge is to bring back the HVDC 
grid, e.g. the DC voltages, to a normal operation state after 
a contingency have happens. 

Stability analysis of combined AC and MTDC systems was 
performed in 0. in m and 0, decentralized controllers 
are employed to share frequency control reserves. In 0 
no stability analysis is performed, whereas IS guarantees 
stability provided that the connected AC areas have identical 
parameters and the voltage dynamics of the HVDC system 
are neglected. EOl considers an optimal decentralized con¬ 
troller for AC grids connected by HVDC systems. 

By connecting the AC areas with a communication net¬ 
work supporting the frequency controllers, the performance 
can be further improved, compared to a decentralized con¬ 
troller structure without such communication. In this paper, 
we seek to explore controllers which improve performance 
of existing controllers. For this, we first propose a controller 
performance measure. 

Several distributed and decentralized controllers for shar¬ 
ing frequency control reserves have been proposed in the 
literature. In HE a distributed controller, relying on a 
communication network, was developed to share frequency 
control reserves of asynchronous AC transmission systems 
connected through an MTDC system. However, the controller 
requires a slack bus to control the DC voltage, and is 
thus only able to share the generation reserves of the non¬ 
slack AC areas. Another distributed controller is proposed in 
m. Stability is guaranteed, and the need for a slack bus 
is eliminated. The voltage dynamics of the MTDC system 
are however neglected. Moreover the implementation of the 
controller requires every controller to access measurements 
of the DC voltages of all MTDC terminals. In fl3l . fT4l 
distributed secondary generation controllers are proposed. 


where the MTDC dynamics are explicitly modeled, and the 
DC voltages are controlled in addition to the frequencies. The 
controller doeas not rely on a slack bus for controlling the DC 
voltages. The distributed control architecture is more scalable 
than a centralized architecture where information from all 
controllers has to be processed simultaneously. By using local 
and neighboring state information, we propose controllers, 
which can be implemented even when communication is 
unavailable. This paper builds on the results in fl5l . 1 1 3j , 
QD, but significantly generalizes the models of the power 
system. 

The remainder of this paper is organized as follows. In 
Section [II] the system model and the control objectives are 
defined. In Section|IIIJ a distributed secondary frequency con¬ 
troller for sharing frequency control and restoration reserves 
is presented, and is shown to satisfy the control objectives. In 
Sections [TV] and [V] the results are generalized to general AC 
networks and 7r-link models of the HVDC lines, respectively. 
In Section [Vi] simulations of the distributed controller on a 
six-terminal MTDC test system are provided, showing the 
effectiveness of the proposed controller. The paper ends with 
concluding remarks in Section [VH| 

II. Model and problem setup 

A. Notation 

Let Q be a static, undirected graph. Denote by V and £ 
the set of vertices and edges of Q , respectively. Let M, be 
the set of neighboring vertices to i £ V. Denote by Cw the 
weighted Laplacian matrix of Q , with edge-weights given by 
the elements of the diagonal matrix W ED- Let e* denote 
the ith Cartesian unit vector. Let C' denote the open left 
half complex plane, and C _ its closure. We denote by c nxm 
an n x to- matrix, whose elements are all equal to c. For 
simplifying notation, we write c n for c nx 1 . 

B. Model and objective 

We will give here a unified model for an MTDC system 
interconnected with several asynchronous AC systems. We 
consider an MTDC transmission system consisting of n 
converters, denoted i = 1 ,n, each connected to an AC 
system, i.e., there are no pure DC nodes of the MTDC grid. 
The converters are assumed to be connected by an MTDC 
transmission grid, i.e. there exist only one connected MTDC 
grid and not several MTDC grids. The node connected to 
converter i is modelled by 

OiVi = - J2 iriVi-vA + iF, (i) 

jeNi "" 

where V, is the DC voltage of converter node i, Ci > 0 
the total capacitance of the converter and the HVDC line 
connected to the considered converter, and I'f the injected 
current from the DC converter to the DC node. The constant 
Rij denotes the resistance of the HVDC transmission line 
connecting the converters i and j. The MTDC transmission 
grid is assumed to be connected. Note that the converter 
model 0 of the MTDC system does not take the dynamics of 
the HVDC lines into account, caused by the inductance and 


capacitance of the lines. In Section [VJ however, we show 
that the model (|T[i can be generalized to a 7r-link model, 
where each HVDC line consists of an arbitrary number of 
resistors, inductors, and capacitors in series. Only HVDC 
nodes which are connected to a converter are considered 
in our model 0- This implies that intermediate nodes are 
not captured by the model. Modelling intermediate nodes 
would result in differential-algebraic equations, resulting in a 
far more complex analysis. While systems with intermediate 
nodes can be transformed into systems without intermediate 
nodes by Kron reduction ED, this is beyond the scope of 
this paper. Each AC system is assumed to consist of a single 
generator which is connected to a DC converter, representing 
an aggregated model of an AC grid. The dynamics of the AC 
system are given by (3D: 

multi = Pf’ 3 + Pf - Pf j , (2) 

where m* > 0 is its moment of inertia. The constant Pf n 
is the generated power, Pf is the power load, and Pf 
is the power injected to the DC system through converter 
i, respectively. While the model |2]) is restricted to single¬ 
generator AC systems, we show in Section[IV]that this model 
can be generalized to a network of arbitrary many generators. 
The control objective can now be stated as follows. 

Objective 1. The frequency deviations are asymptotically 
equal to zero, i.e., 

lim uii(t) — w ref = 0 i = 1,..., n, (3) 

£—>•00 

where w ref is the nominal frequency. The total quadratic cost 
of the power generation is minimized asymptotically, i.e., 

lim^oo Pf™ = Pf™*, V* = 1,..., n, where 

[Pr,...,Pri =argmin^/f (Pf n ) (4) 

subject to ([3]), i.e., Pf™ + Pf — Pf =0, Vi = 1,..., n 
and l P," 11 = 0, i- e -, power balance both in the AC grids 
and in the MTDC grid. The positive constants ff represent 
the local cost of generating power. Finally, the DC voltages 
are such that the a quadratic cost function of the voltage 
deviations is minimized asymptotically, i.e., lim^oo Vi = 
V*, Vi = 1,... ,n, where 

1 " 

Wi > • • • , V n] = argmin - ]T ff (V - Vff 2 (5) 

subject to ([3])-([4]), and where the f] is a positive constant 
reflecting the local cost of DC voltage deviations and Vf f is 
the nominal DC voltage of converter i. 

Remark 1. Note that the order in which the optimization 
problems 0-0 are solved is crucial, as 0 and the optimal 
solution of 0 are constraints of 0. 

Remark 2. The minimization of ([4]) is equivalent to power 
sharing, where the generated power of AC area i is asymp¬ 
totically inverse proportional to the cost ff. The cost ff 


can be chosen to reflect the available generation capacity of 
area i. 

Remark 3. It is in general not possible that lim^oo V. (t) = 

ref 

V i 1 Vi = 1,..., n, since this does not allow for the currents 
between the HVDC converters to change by 0 - Note that the 
optimal solution to (|4| fixes the relative DC voltages, leaving 
only the ground voltage as a decision variable of < 0 >- Note 
also that the reference DC voltages Vf e ^, i = 1 ,n, are 
generally not uniform, as is the reference frequency U) re f 

Remark 4. Note that Objective [ 7 ] does not include con¬ 
straints of, e.g., generation and line capacities. This requires 
that the perturbations from the operating point are sufficiently 
small, to guarantee that these constraints are not violated. 

III. Distributed frequency control 
A. Controller structure 

In this section we propose a distributed secondary fre¬ 
quency controller. In addition to the generation controller 
proposed in EL we also propose a secondary controller 
for the power injections into the HVDC grid. We implement 
the controllers for single AC generators. In Section [TV] we 
generalize the controller for AC grids of arbitrary size. 

The distributed generation controller of the AC systems is 
given by 

7 SV 

PT = -K?°° v {ui - ^ ref ) - ^f oop - \ 

Vi = K? oop '\u]i - w ref ) - - rjj), (6) 

jeNi 

where Kf°° v and Kf °° v ‘ 1 are positive controller parameters. 
Moreover, c r h = cl > 0, i.e., the communication graph 
is supposed to be undirected. The above controller can be 
interpreted as a distributed Pi-controller, with a distributed 
averaging filter acting on the integral states //,. The first line 
of Equation ([6]) resembles a decentralized droop controller 
with a setpoint given by rp. The second line of Equa¬ 
tion ([6]) updates the variable rp in a distributed fashion by 
a distributed averaging integral controller. The magnitudes 
of the variables c7 determine how fast the generated power 
levels converge. While a larger magnitude of c7- could lead 
to faster convergence of the generated power, it can also 
induce oscillations. It is possible to implement a decentralized 
version of ([6]) by dropping the states rp . This results in the 
following controller 

Pf n = -/Tf°°P( Wi - w ref ). (7) 

The proposed converter controllers governing the power 
injections from the AC systems into the HVDC grid are given 
by 

if J = K?( Ui - w ref ) + KY{Vf* - V.) 

+ 4(0* - 

j&Ni 

K UJ 

0* = Wi ~ 70*, (8) 


where K] and Kf are positive controller parameters, and 
if’ nom is the nominal injected power, 7 > 0 and f = f > 
0. If 7 = 0, the converter controller ([ 8 } can be interpreted 
as an emulation of an AC network between the isolated 
AC areas, as it resembles the swing equation. The auxiliary 
controller variables <fii are then equivalent to the phase angles 
of AC area i, whose differences govern the power transfer 
between the areas. Larger magnitudes of cf correspond 
to higher conductances of the AC lines, and thus stronger 
coupling and faster synchronization of the frequencies. If 
7 > 0 , damping is added to the dynamics of fi. Damping 
generally improves stability margins, and turns out to be very 
useful in the stability analysis. However, a nonzero 7 also 
implies that the AC dynamics are not emulated perfectly. This 
implies that exact frequency synchronization might not be 
possible in general. In contrast to a connection with AC lines, 
the power is fed into the MTDC grid and then transfered 
to the other AC areas through the MTDC grid rather than 
through an AC grid. Also the converter controller can be 
implemented in a decentralized version by dropping the states 
<j>i, resulting in the following controller 

Pf = Kf ( Ui - u/ ef ) + KY {Vf ei -Vi). (9) 

The HVDC converter response is assumed to be instanta¬ 
neous, i.e., injected power on the AC side is immediately 
and losslessly converted to DC power. This assumption is 
reasonable due to the dynamics of the converter typically 
being orders of magnitudes faster than the primary frequency 
control dynamics of the AC system J4]]. The relation between 
the injected HVDC current and the injected AC power is thus 
given by 

Vjf ] = if. (10) 

By assuming Vj = V nom i = 1,..., n, where V Dom is a global 
nominal DC voltage, we obtain 

ynomjinj _ pin] ( 11 ) 

Assumption GD relies on the assumption that the voltages V t 
do not deviate significantly from the nominal voltage V nom . 
Since for most HVDC converters the acceptable deviation 
from the nominal voltage is less than 5% ED, the approxi¬ 
mation mi would result in a relative error smaller than 5%. 

B. Stability analysis 

We now analyze the stability of the closed-loop system. 
Define the state vectors Cj = oj — w ref l„ and V = V — 
V iei , where w = [w 1 , - - -, w„] T , V = [Vi,, V n ] T , V ief = 
[V{ et ,..., Vf et ] T , V = [Vi, • • •, Vn] T , and 0 =[</>!,..., 0„]. 
Combining the MTDC {T}, the AC dynamics ([2]) with the 
generation control ([ 6 ]), the converter controller 0 with the 
power-current relationship (fTT), we obtain the closed-loop 


dynamics 

Cj =m( - (iv droo P + k u )6j + K v V 

—K v (A ' w )" 1 A' droop ’ \ - C^fi + P m ) 

x 1 ( K v \ - 1 

H=-—AAT“w - E \C R + — — V + -— EC^ 

T/nom l - rL T/nom / T/nom 

7? =A -droop,I ( 2 I _ 

fi ={K V )- X K U (C - -70, ( 12 ) 


where M = diag(mi _1 ,..., m„ _1 ) is a matrix of inverse 
generator inertia, E = diag(C 1 " ,..., Cf 1 ) is a matrix of 
electrical elastances, C R is the weighted Laplacian matrix 
of the MTDC grid with edge-weights 1 / R l3 , C v and C q j 
are the weighted Laplacian matrices of the communication 
graphs with edge-weights </ and cfj, respectively, and 
pm _ jprn^ _ _, P™] T - We define the diagonal matrices of 
the controller gains as AT W = diag(iL^,..., Kf), etc. 

Let y = [lu t ,V t ]' 1 define the output of ( fl2| ). Clearly the 
linear combination 1 P <p is unobservable and marginally stable 
with respect to the dynamics ( fl 2 l >, as it lies in the nullspace 
of €(/,. In order to facilitate the stability analysis, we will 
perform a state-transformation to this unobservable mode. 
Consider the following state-transformation: 




where S is an n x (n — 1 ) matrix such that 4= l n 
orthonormal. By applying the state-transformation © to 
©, we obtain dynamics where it can be shown that the state 
(f)\ is unobservable with respect to the defined output. Hence, 
omitting <j>\ does not affect the output dynamics. Thus, we 
define <j>" = [<p' 2> , </>' n \, and obtain the dynamics 


Cj =m( - (A' droop + K u )Cj + I< v V 

—A' y (A'“) “ 1 A' droop ' l rj - C^Sfi" + P m ) 

x 1 ( K v \ - 1 

V = -— EK u Cj -E[C r + -— V + -— EC^Sfi" 

T/nom \ ' T/nom / T/nom ™ p 

f) =K dm °p-Y - c ri n 

j)"=S T (K v )- 1 K u Cj-~j<j>". ( 14 ) 


We are now ready to show the main stability result of this 
section. The following assumptions are later used as sufficient 
conditions for closed-loop stability. 


Assumption 1. The Laplacian matrix satisfies C n , = k rll C R - 

Assumption IT] can be interpreted as the emulated AC 
dynamics of ([ 8 ]» having the same susceptance ratios as the 
conductance ratios of the HVDC lines. Assumption |T] can 
always be satisfied by appropriate choices of the constants 
Cij in ([Sj. 

Assumption 2. The gain 7 satisfies 7 > k t / > /(4y nom ). 

Assumption [2] lower bounds for the damping coefficient 
of the converter controllers. Note that the bound on 7 is 


independent of the topology of the communication network. 
This is particularly desirable in a plug-and-play setting, where 
new nodes can be added to the system, without having to 
change 7 . 

Theorem 1. If A ss umptions [7] and [2] hold, the equilibrium of 
© is globally asymptotically stable. 

Proof: The proof follows from Theorem [4] and is thus 
omitted. ■ 

Corollary 2. Let Assumption U\ hold and let 7 , fc e , be 

given such that Assumption lT holds. Let K v , K w and 
K droop.I be such that (P/-1 = K V( K u)-l K droop,I and 

F v = K v , where F p = diag(/j P ,..., f p ) and F v = 
diagj// ..., //). Then the dynamics © satisfy Objec- 
tiveulin the limit when || (K U )~ 1 K^ / 11^ —>• 0, provided that 
the disturbance Pf l is constant. 

Proof: By Theorem 0 © has a unique and stable 
equilibrium. Letting (ft" = 0„_i implies S T (K v )~ 1 K u co — 
yfi" = 0„_i. Now ||(AT u; )- 1 Ar y || oo —)• 0 implies that 
S t Cj = 0 <£> w = k\\ n for some k\ £ M. Letting 77 = 0„ in 
© yields 

A' dr °oP- - C v r, = 0„. 

By inserting Cj = kil n and premultiplying the above equation 
with 1/ we obtain that k\ = 0, so Co = 0 n so Equation ([3]) 
of Objective 0 is thus satisfied. Thus 77 = k 2 \ n for some 
k 2 £ R. Finally, we let V = 0 n in ( fl4| >: 

K u Cj - (v nom C R + K v \ V + C^Sfi" = 0 n . (15) 

Inserting Cj = 0„ and premultiplying ( p~5j ) with / yield 

l p K v V = 0. (16) 

Inserting Cj = 0„ and 77 = k 2 l n in *0 yields 

P gen = -k 2 K v (K u )~ 1 K dTOOp ’ I l„, (17) 

where P gen = [Pf™, • ■ ■, P ? en ] 2 . It now remains to show that 
the equilibrium of ( |T4l > minimizes the cost functions <[4]) and 
0 of Objective 0 Consider first 0, with the constraints 
Pf n + - P/ = 0 , i = 1 ,..., n and /” =1 P“ j = 0. 

By summing the first constraints we obtain P® en = 

- EILi ^T- The KKT condition of 0 is 

P P P gen = -fc 3 l n . (18) 

Since (P p ) _1 = K v (K^Y ' A' dr00p ’ ! , ([T7| and ([18]) 

are identical for k 2 = Y. We conclude that 0 is 
minimized. Since P gen = -K V (K U )- 1 AT droop ’ = 

—k 2 K v (K^Y 1 K dmop - l l n and Cj = 0 n , premultiplying the 
first n rows of the equilibrium of ( p~4] ) with M -1 , and adding 
to the (?i + l)th to 2 ?ith rows premultiplied with )/ nom A _1 
yields 

-F nom C R V - k 2 K v {K u )" 1 AT droop ’ 1 1 „ = P m . 

Premultiplying the above equation with l p yields k 2 = 

- E"=i P r E"=i k vY<°°v ■ Additionally, C R V is uniquely 















determined. Now consider ([5]). Note that P" 1J and hence J‘ nj , 
are uniquely determined by 0. By the equilibrium of Q, 
C R V = I m J, where P n « = [Iff ..., I" ,J ] T . Thus, the KKT 
condition of ([5]) is 

F v V = C R r, (19) 

where r £ R". Since C R V is uniquely determined, we 
premultiply ( fl9| with 1 f and obtain the equivalent condition 

1 T n F v V = 0. (20) 


A. Objective 

Consider the AC transmission grid connected to converter 
i, and suppose it consists of m generator buses. Without 
loss of generality, we may assume that converter i of the 
MTDC grid is connected to generator i\ of the AC system 
i. Let d, k be the phase angle of bus i\.. The dynamics of 
the power system are assumed to be given by the linearized 
swing equation ED, where the voltages are assumed to be 
constant. As before, we consider the incremental states with 
respect to their reference values: 


Since F v = K v , ( fT 6 | ) and ( |20| are equivalent. Hence ([5j is 
minimized, so Objective [I] is satisfied. ■ 

Remark 5. Corollary [2] provides insight in choosing the 
controller gains of 0 and 0. to satisfy Objective |7] 

While the generation controller 0 and the converter 
controller 0 offer good performance in terms of satisfying 
Objective [T] it may not be possible to implement these 
distributed controllers, e.g., due to lack of communication 
infrastructure. For such MTDC systems where appropriate 
communication is lacking, it may be desirable to instead 
implement decentralized generation and converter controllers. 
In other situations it might be possible to implement the dis¬ 
tributed generation controller 0 , while it is more desirable 
to have the HVDC converters operating independently with 
decentralized controllers. In the following corollary, we show 
that the decentralized generation and converter controllers 0 
and 0 also globally asymptotically stabilize the combined 
MTDC and AC system. 

Corollary 3. Let Assumption |T] hold and let 7, k r ., be given 
such that Assumption [2] holds. Consider the dynamics of the 
MTDC dynamics 0 and the AC dynamics 0 with the gen¬ 
eration controller 0 or 0 , respectively, and the converter 
controller 0 . The equilibria of the resulting closed-loop 
systems are globally asymptotically stable. 

Proof: The proof is in line with the proof of Theorem [4] 
where we discard the variables 77 and <fi. ■ 

While the optimality results of Corollary [2] do not hold 
for any other controller combinations than 0 and 0 , the 
following remark can be made about the average frequency 
errors. 


Lemma 1. Consider the dynamics of the MTDC dynamics 
(1 1 and the AC dynamics 0 with the generation controller 
(61 and the converter controller 0 . Any equilibrium of 
the resulting closed-loop system satisfies Ff oop ' ! { oj 1 — 

uj re f) : 0 , i.e., the average frequency errors are zero. 


Proof: Consider the closed-loop dynamics ( fl2| ). Letting 
77 = 0 „ and premultiplying this equation with 1 J yields 0 „ = 

l^oop-'w - 1^77 = £” =1 Af oopJ (wi - w ref )- ■ 


hi. k — ! -C'i 


m ik ^ ik =-{Kfr + Kl)u: ik - 




'y ] hikj(hik 




( 21 ) 


where S lk is the phase angle and Chi k = uji k — u/ ef is the 
incremental frequency at bus ik, rri lk > 0 is the inertia 
of bus ik, ki k j = where Vi is the constant 

voltage of bus i, and bi k j is the susceptance of the power 
line Moreover A'R°° P = 0 for k ^ 1 , since power 

injection through the HVDC converter only takes place at 
bus i\. The constant Pf cn is the generated power by the 
generation control, Pjj is the uncontrolled deviation from 
the nominal generated power at generator ik, respectively. 
The variable P™ J = 0 for k f 1 is the power injected to 
the DC system through converter i. We assume that the AC 
voltages are constant, thus implying that fcy is constant. In 
order to account for the additional generators, we need to 
slightly modify Objective [T] 

Objective 2. The frequency deviations converge to zero, i.e., 


lim uji k (t) — u/ ef = 0 k = 1 ,... 7 ij, i = 1 ,..., n. ( 22 ) 


The total cost of the power generation is minimized asymp¬ 
totically, i.e., linii^oo Pf en = P ® en *,i = 1,... ,n, where 


[prv-.jr*] 


1 

= argmin - 

Pi P n ^ 




(23) 


subject to 1 ^. (P® en + P™ — P“- i ) = 0,i = 1 , ...,n and 
2"=i ^ mJ = 0, i.e., power balance both in the AC grids 
and in the MTDC grid. Here Jf en = [Pf™,..., Pf“] T , i = 
1,... ,n. Finally, the DC voltages are such that a quadratic 
cost function of the voltage deviations is minimized asymp¬ 
totically, i.e., lim^oo Vi = V*, i = 1 ,..., n, where 

n 

[Hr,.. •, v:] = argmin - £ fY(V t - V^) 2 (24) 

subject to 0 - 0 . Here ff and fY are positive constants. 


IV. Generalisation to AC generation network b. Controller structure 


In this section we generalize the single-generator model of 
Section II-B to an AC grid with arbitrary size. 


In this section we generalize the distributed secondary 
frequency controller 0 and the converter controller 0 to 





the full AC network. The distributed generation controllers 
of the AC network i are in this case given by 


K) 


pgen _ _ T^droop- _ 




Vi = Y K ik° V,lCj i k - Y c ij (Vi -Vj)- (25) 

jeM'i 


fe =1 


where A'‘*‘ oop , Kfa, Kf and A l droop ’ 1 are positive controller 
parameters, and c0- = c' , ; > 0. Compare Equation <[ 6 ]>. 
The above controller can be interpreted as a distributed PI- 
controller, with a distributed consensus filter acting on the 
integral states 77 . The converter controller governing the 
power injections from bus i\ of the AC system i into the 
HVDC grid is given by 

pinj = pint, nom + ^ _ ^ref) + R V fyv* _ y.) 

+ Y c tfafa - fa) 

i&Mi 

fa = “ wref ) “ 70i, (26) 

where 7 > 0 and cfj = <fa % > 0. Compare Equation ([ 8 ]). In 
vector-form pT) becomes 

Si — UJi 

fa = Mfa — (Af oop + K?)£* - Cf8i 

+ Pf en + P™ — Pfafa , (27) 

where S t = [6 h ,..., S irii } T , uj z = [c u^,... ,u in .] T , 

Mi = diag Cf c is the Laplacian ma¬ 
trix of the graph corresponding to the AC transmis¬ 
sion system, with edge-weights given by ki k j, Kf loop = 
diag(< oop ,...,<° op ), K? = diag(K%,0,...,0), and 
Pf cn = [Pffa ,..., Pf en ] T , etc. Consider the output yi = (bi of 
( |27) . With respect to y, = uj, , the dynamics have a marginally 
stable unobservable mode. Thus, similar to Section |III-B| we 
consider the state transformation 






1 1 

E*h e 

1_1 

w; n 

Si 

fa 

fa = 


where Si is an i n x (i n — l)-matrix such that [4=1„, Si] is 
orthonormal. It can be shown that V is unobservable, and 

i l 

can be omitted by introducing the state 5" = [61,... ,8fa ] 1 . 
This state-transformation results in the dynamics 

j" = S?u>i 

fa = Mi ( - (ATf oop + Ki)Cji - CfSiS'f 

+ pge n + pm_pmj)_ (28) 

Since the input-output dynamics of ( |27| and ( |28) are iden¬ 
tical, we henceforth only consider the dynamics ( |28) . By 
combining the dynamics Q and ( |28T l with the controllers 
(|25) and ([26]), and considering the change of coordinates (p~3) 


and <j>" = \(j>fa ..., (f>' n ] we obtain the dynamics 
5" =Sj'u>i, i = 1 ,...,n 

fa =Mi ( - (Af oop + K?)Qi + e-jfay fa - CfStf 

'I niVi - er + P™), i = 1 ,... ,n 
A 1 / K v \ - 1 

V = -—EK uj Cj -E\C r + —— V + — — EC^Sfi" 

T/nom l 1 T/nom / T/nom p p 

n 

V =Y eil n K i tOOP ’ l fa- C vV 

i =1 

fa'=S T (K v )- 1 K u ui - 74-10", (29) 

where u = ,..., w ni ] T , ATf oop = 

diag(A^ roop ,...,A^° op ), Af oop ' 1 

diag(<°° p - ht °° P ’'). K“ = diag ,..., Kffa. 

Theorem 4. The equilibrium of the dynamics ( |29[ > is globally 
asymptotically stable under Assumptions [7] and [2] 

Corollary 5. Let Assumption [7] hold and let 7 , k n , be given 
such that Assumption [ 2 ] holds. Let I\fa. Kf and Kfa oop be 
such that (Ff)- 1 = Kfa{Kf)- 1 Kf' oop , i = 1 ,...,n 
and F l = K l , where F p = diag(/4 • • •, f, f) and 
F p = diag(/ ] l/ ,..., ffa ). Then Objective [2| is satisfied in 
the limit when || (K U1 )~ 1 K V || —> 0, provided that the 
disturbance P™ is constant. 

Proof: Consider ( |29| . Letting 5" = 0 ni i = 1 
yields u), = ki\ Ui i = 1,..., n. Letting <j>" = 0 n yields 

S t (K v )- 1 K“u> - 7 J n _r = 0„. 


Now ||(A' W )- 1 AT V || —> 0 in the above equation implies 
S t Cj = 0 G) = kl n , k £ K. This implies u>i = kl ni ,i = 
1 Letting r) = 0 n , inserting £7 = kl n . i = 1 ,n, 

and premultiplying the equation with l p finally yields A; = 0, 
and thus £7 = 0 ni ,i = 1,... ,n, i.e., ( |22] i is satisfied. Letting 
fj = 0 n and inserting i2>i = 0 ni i = 1 ,..., n yields rj = kil n , 
which inserted in ([25) yields 


K v 

pS en _ l. i p-droop, I 
1 TSUJ 

h 


k = 1 ,..., rii- 


(30) 


Linally we let V = 0 n , insert lu = (),, and premultiply the 
equation with l^C and obtain 

l p K v V = 0. (31) 

By similar arguments as in the proof of Corollary [2] we can 
show that ( |3()) and ( |3T) are equivalent to the KKT conditions 
of (|23) and (|24), respectively. This concludes the proof. ■ 


V. Generalization to 7t-Link HVDC model 

In this section we extend the HVDC line model to consider 
the inductance and capacitance of the HVDC lines. We model 
the HVDC lines as series of £ 7 r-links consisting of resistors, 











inductors, and capacitors. The dynamics of HVDC line k, 
connecting converters i and j, are given by 

CiVi = C - 51 4,i 

fcejVf 

-4-4, 1 = —-Rfc4,l + Vi — Vk,i 
Cl m \ 4,r = 4,i - 4,2 
LkR,2 = ~ 44,2 + 4,1 — ^4,2 


4fe ne Vfe,f-i = 4,1-t — h,t 

Lklk,t = —Rkh,e + Vk,e- 1 — Vj 

** - 1 - E 'w. 

fee _V° U ' 


/; 


(32) 


where 14, g and 4,g denote the DC voltage and current 
of line segment <y, respectively, and G'j] nc , /■?/,,, and 4 
are the capacitance, resistance and inductance of each line 
segment of line k, respectively. The sets Aff and J\f° ut 
denote the incoming and outgoing HVDC lines to converter i, 
respectively. To simplify the derivations we only consider AC 
areas consisting of single generators, but the results can be 
generalized to also include the AC generator network model 
of Section IV Combining the DC voltage dynamics (J32]», the 


frequency dynamics |2]i with the generation control (|6|, the 
converter controller (|8]» with the power-current relationship 
HD. with <j>" debited as in Section [TTTI we obtain the closed- 
loop dynamics 


&=M( - (/C droop + K u )Cj + K v V 

—K v (K u )~ 1 K droop ’ 4? - C^Sf" + P m ) 

(~ KV V ~ ^ n ° m An4 + V n ° m D out I £ 
+K“u> + C^Sf") 
ii=L~ 1 (—RI 1 + D?V — Vi) 

Vi=P line (4 - 4) 

i 2 =L-\-RI 2 + V 1 -V 2 ) 


I4_i=P line (4-i - 4) 
i^L-'i-Rh + Vi-t- D oui V) 
ri=K imop - 1 uj - £ v r) 

<j>"=S T (K v )- 1 K“ut- 74-10", (33) 

where M, E K dmop , K u , K v , K dmo p- t , Cs, C r , 

P m , ih, V are debned as in Section |Hl| P llne = 
diag((Ci me ) _1 ,... , (G^ ne ) _1 ), L = diag(Li,.. .,L m ), and 
R = diag(f?i,... ,4 m ). The n x m-matrices D; n and D out 
describe the adjacency relations between the in- and outgoing 
HVDC lines and the converters. Element (i,j) of l) m or P out 
is 1 if line j originates or terminates at converter i, and 0 
otherwise. The variables I q = [4, 9 — 4° m ,..., I m , q — J™ m ] T 
and V q = [14 , q - Vf° m ,..., V m , q - V™”] T debne the 
incremental currents and voltages of the q:th line segment 


for all HVDC lines 1 ,,m. The constant /" onl denotes the 
nominal current of each line segment of line k, and Vj?° m 
denotes the nominal voltage of segment q of line k. 

Theorem 6. The equilibrium of the dynamics HD is globally 
asymptotically stable under Assumptions [7] and [2] 

VI. Simulations 

In this section, simulations are conducted on a test system 
to validate the performance of the proposed controllers. The 
simulation was performed in Matlab, using a dynamic phasor 
approach based on ED- The test system is illustrated in 
Fig-0 We model the HVDC line as a single 7r-link, with 



Fig. 1: Test grid consisting of 6 asynchronous IEEE 14 AC grids connected 
through an MTDC grid. 


parameters given in Table [I] The terminal capacitances are 
Ci = 0.375 x 10“ :1 p.u. The AC grid parameters were 
obtained from ED. The generators are modeled as a 6th 
order synchronous machine model controlled by an automatic 
voltage controller and a governor 01- The consumers in all 
AC grids are modeled as constant power loads. It is assumed 
the linear power-current relation (flOf holds. 


Table I: HVDC grid line parameters 


(hi) 

Rij [p.u.] 

L i:j [10 3 p.u.] 

Cij [p.u.] 

(1, 2), (1, 3), (2.4). (3,4) 

0.0586 

0.2560 

0.0085 

(2, 3) 

0.0878 

0.3840 

0.0127 

(2, 5), (4,5) 

0.0732 

0.3200 

0.0106 

(2, 6), (3, 5), (5. 6) 

0.1464 

0.6400 

0.0212 




Table H: 

Controller parameters 


Kf 

KV 

7 ^-droop 

Tv-droop, I 

4 

C - ■ 

1-3 

4> 

C- *7 

ij i 

1501 

80 

9 

3.35 

5 /Rij 

15 /Rij 0 


































































The three different controllers proposed in this paper 
are applied to the test grid, i.e., © and 0 . 0 and 0 , 
and 0 and 0, with parameters given in Table [II] The 
communication network of 0 and 0 is illustrated by the 
dashed lines in Fig. [T] 

We set 7 = 0, so Theorem [4] does not guarantee stability 
of the equilibrium. However, the closed-loop system matrix 
can easily be verified to be Hurwitz. At time t = 1 the 
output of generator 2 in AC area 1 was reduced by 0.2 p.u. 
Fig. [2] shows the average frequencies of the AC grids for all 
three controllers. Immediately after the fault the frequencies 
of the AC area of the increased load drop. The frequency 
drop is followed by a DC voltage drop in all converter 
nodes, due to © or 0 . We note that the frequencies are 
restored to the nominal frequency for © and © as predicted 
by ( |22| >. However, for the controller combination 0 and 
© the frequencies are not restored to the nominal values. 



Fig. 2: Average AC area frequencies. After a reduced generation of 0.2 
p.u. at t = Is, the frequencies synchronize fast, and are subsequently 
restored to the nominal frequency when the distributed generator and 
converter controllers {6} and j8]i are employed. For the controllers JbJ and 
0 or 0 and 0. a static control error remains present. 



Fig. 3: DC converter voltages. After a reduced generation of 0.2 p.u. 
at t = 1 s, an immediate voltage drop is followed by a restoration of 
the average voltage errors to zero when the distributed generator controller 
(6} is employed together with either converter controller or |9}. This 
corresponds to the minimization of {5}. However, when is employed, the 
average voltage error does not converge to zero. 



Fig. 4: Total increase of generated power in the AC areas. After a reduced 
generation of 0.2 p.u. at t = 1 s. The generated power is increased for 
all AC areas regardless of the controllers employed. When the distributed 
generator and converter controllers {6} and 0 are employed, all generators 
generate to the same power and have the same marginal generation costs 
asymptotically, thus minimizing {4} asymptotically. For 0 and 0 or 0 
and 0, the marginal generation costs do not converge to the same value. 


despite that a secondary frequency controller is employed. An 
intuitive explanation to this is that the distributed frequency 
controller © requires that the frequencies in the different 
areas synchronize. This synchronization is achieved by ©, 
but not by ©. 

Fig. a shows the DC voltages of the terminals. A similar 
behavior as from the frequencies can be seen. As foreseen, 
the voltage deviations can not be controlled back to zero, 
since this would result in zero current flows. For the con¬ 
trollers 0 and ©, the average voltage deviation is restored 
to zero as predicted by ( |Ti~| ). 

Fig. [4] shows the total increase of the generated power 
within each AC area. Initially all generators have a similar 
oscillating behavior, due to the AC voltage oscillations. For 
the controllers © and © the total increase in the generated 
power of the AC areas converge to the same value, thus 
minimizing ( |23) l, since the controller parameters are chosen 
uniformly for all converters i. For the other two controller 
combinations, the increase in the generated power of the AC 
areas does not converge to the same value. This implies that 
( |23| ) is not minimized, and we conclude that the simulations 
confirm that the controllers 0 and © have superior power 
sharing properties. Clearly, the communication in the con¬ 
trollers 0 and © is essential to eliminate the frequency 
error and sharing frequency control reserves. We can note 
relatively high frequency oscillations in the transient response 
of the generated power. These oscillations originate from the 
averaging term in the dynamics of the internal controller 
variables rji. By reducing the constants cT, these oscillations 
can be reduced at the expense of slower convergence. 

Fig. 0 shows the power set-points of the converters. For 
controller 0 and © the control in area 1 overcompensate 
the disturbance at the beginning, also the other areas have 
some minor oscillations. A stable operation point is found 
after 6 s. Area 2-6 have all an equal share of the disturbance. 
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Fig. 5: Converter power set-points for all 6 areas. After a reduced 
generation of 0.2 p.u. at t = 1 s. For all controllers the area with the 
disturbance starts to extract power from the MTDC grid and all others inject 
power. When the distributed generator and converter controllers JbJ and {8} 
are employed, the converter set-points stabilizing after 6 s to a point with 
perfect sharing of the disturbance. For 0 and j9j or {7} and it takes 
some more time to find a stable point and the converters differ from each 
other. 


It takes longer time for the other controllers to find a stable 
operation point, where controllers ([ 6 ]) and (JyJi are significantly 
faster than |7]) and (|9j. The sharing of the disturbance are for 
both controller not equal, indicated by the different power 
set-points of the converters. 

VII. Conclusions 


Differentiating W along trajectories of ([29]) yields 


n ffu t 

W = E W WfSjCfSiSjCH 

i—1 i \ 

.-droop 


+ ( - (Af oop + Kf)CH + er KY V t - CfSrf'f 

-^K^’ l l ni r li -e 1 eJC,Srfj 

+ V T - ( V nom C R + K v ) V + C+Stf'') 

+ e i l n K ?°° V ’ l(b i - 4 ^ 

+ <t>" T S T /: 4l s(s T (K v )- 1 K“u - 7 / n _ 1 0") 

Since SS T = I n - 4 and S z Sf = I n . - ^1 we 
have C^SS T = C$ and Cf c SiSf = Cf c ll22ll . Furthermore 

11 TSUI 

Y eief =u T (K v )~ l C^Sf" 


i—1 1 


A'j 1 


E KYVi=u T K“V 

i—1 i 


V T Y ^nXf^^Y ^ X n K \ 


T 7^droop,I 




i—1 


By defining 

V' = 


-4 In S 


i—1 


V V = 


-4 In S 


v\ 


In this paper we have studied distributed secondary con¬ 
trollers for sharing frequency control reserves of asyn¬ 
chronous AC systems connected through an MTDC sys¬ 
tem. The proposed controllers were shown to stabilize the 
interconnected AC systems and the MTDC grid. The AC 
grid frequencies were shown to converge to the nominal fre¬ 
quency. Furthermore, quadratic cost functions of the voltage 
deviations of the MTDC terminals and of the generated power 
was minimized asymptotically. The results were first derived 
for single-generator AC grids and purely resistive MTDC 
grids, and later generalized to AC grids of arbitrary size, and 
HVDC links modelled as 7 r-links. Finally, the results were 
validated on a six-terminal MTDC system with connected 
IEEE 14 bus AC grids. Future work will focus on stability 
results under communication delays. 


Appendix 


Proof of Theorem [7J Without loss of generality, let 
If" = 0 n ., i = 1 ,,n. Consider the Lyapunov function 
candidate, which is positive definite and radially unbounded: 


n k u / \ 

W-- YjfY ( W) T S?C?Stf + 

V n ° m ~ 1 rj, 1 llT T 

+ —V T CV + -ifr, + -<j>" T S T C 4 ,S<t>"- 


we obtain VjC R V = V" T S T C R SV" and = 

V" T S T CtjyScf), where V" = [V.f ..., Vf\ T . We thus obtain 


W < - 


r -| 

T r 

UJ 


V 


L J 



-K u 


-k u 


ib 

K v 


V 



L - 1 


=Q i 


V" 

T 

V nom S T £ R S 

~^S t C r S 

V" 

r. 


- k ^S T C R S 

7 k <t> S T £ R S_ 

r. 


—v'— 

= Q2 


1 

2 


Y^Kt IOOP ^-V T C vV 

i= 1 


where AT droop = diag(A'* oop ,..., 4^ r ° op ). By applying the 
Schur complement condition for positive definiteness, we 
see that Qi is positive definite, since K U (K V )~ 1 (K U1 + 
A'droop) _ AT <j (AV) _ 4\“ = K UJ (K V )~ 1 A' droop > 0. 
By similar arguments Q 2 is positive definite iff ( 7 k^ — 
k‘^/(4V n ° m ))S 1 C R S > 0. Clearly the above matrix inequal¬ 
ity holds under Assumption i since S T C R S > 0, and 
Sx 7 ^ kl n for k 7 ^ 0. Thus W < 0 under Assumption [ 2 ] and 
the set where W is non-decreasing is given by G = {<5" £ 
R ni , i = 1, ..., 71,77 = kl n }, for any k £ R. The largest 
invariant set in G with respect to ( |29| is the origin, since 
Cf c SiS” + ^Kf°° p l l ni rii = 0 Ui implies 8" = 0 Ui and 
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rji = 0. By LaSalle’s theorem, the origin of ( |29l > is globally 
asymptotically stable. ■ 

Proof sketch of Theorem [6[ Without loss of generality, 
let P m = 0,, in ( |33 [ > . One can verify that 

T/nom / ^ t—1 \ 

W=—~ V T CV+ Y Y VfC^Vi 

\ 2—1 2=1 J 

+ 1 -Cu t K“(K v )- 1 M~ 1 Q + * fj T fj + ^ T S T C^Sf 

is a Lyapunov function for m where C lme = 
diag(C'} lne ..., C^ e ). Differentiating W along trajectories of 
( [33] ) where P m = 0 n yields W < 0. By LaSalle’s theorem, 
the equilibrium of ( |33| ) is globally asymptotically stable. ■ 
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